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On Type la Supernovae From The Collisions of Two White 
Dwarfs 
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ABSTRACT 

We explore collisions between two white dwarfs as a pathway for making Type la 
Supernovae (SNIa). White dwarf number densities in globular clusters allow 10 — 100 
redshift < 1 collisions per year, and observations by (Chomiuk et al. 2008) of globular 
clusters in the nearby SO galaxy NGC 7457 have detected what is likely to be a SNIa 
remnant. We carry out simulations of the collision between two 0.6 M white dwarfs 
at various impact parameters and mass resolutions. For impact parameters less than 
half the radius of the white dwarf, we find such collisions produce « 0.4 M of 56 Ni, 
making such events potential candidates for underluminous SNIa or a new class of 
transients between Novae and SNIa. 

Key words: hydrodynamics - nuclear reactions, nucleosynthesis, abundances - 
(stars:) white dwarfs - (stars:) supernovae: general. 



1 INTRODUCTION 

Type la supernovae (henceforth SNIa) play a key role in 
astrophysics as premier distance indicators for cosmology 
(Phillips 1993; Riess et al. 1998; Perlmutter et al. 1999), as 
direct probes of low-mass star formation rates at cosmologi- 
cal distances (Scannapieco et al. 2005; Mannucci et al. 2006; 
Maoz 2008) and as significant contributors to iron-group 
elements in the cosmos (Wheeler et al. 1989; Timmes et 
al 1995; Feltzing et al 2001; Strigari 2006). Our current un- 
derstanding is that there are two major progenitor systems 
for these events. The first possibility, the single-degenerate 
scenario, consists of a carbon-oxygen white dwarf in a bi- 
nary system evolving to the stage of central ignition by mass 
overflow from a low-mass stellar companion (Whelan & Iben 
1973; Nomoto 1982; Hillebrandt k Niemeyer 2000). The sec- 
ond possibility, the double-degenerate scenario, consists of 
the merger of two white dwarfs in a binary system (Iben 
& Tutukov 1984; Webbink 1984; Yoon et al 2007). It is un- 
known at what relative frequency both of these channels 
operate (Livio 2000; Maoz 2008). 

Collisions between two white dwarfs, are likely to hap- 
pen less frequently than binary mergers. However, as dis- 
cussed in Timmes (2009) and Rosswog et al (2009), they 
will occur in globular clusters where the stellar densities are 
extremely high. For a typical globular cluster velocity dis- 
persion of w 5-10 km s _1 , and a white dwarf escape velocity 
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of w 4000 km s _1 , the physical cross-section of equal mass 
white dwarfs is enhanced by a factor of 4 x 10 4 — 2 x 10 5 by 
gravitational focusing. For an average globular cluster mass 
of 10 6 M (Brodie & Strader 2006), a Salpeter IMF (Salpeter 
1955), and a globular cluster core radius of 1.5 pc (Peterson 
& King 1975), we conservatively estimate an average white 
dwarf number density in globulars to be « 10 4 pc -3 , and 
estimating that about 1/1000 of the average 5 x 10 8 M© 
per Mpc -3 stellar mass density in the universe is in glob- 
ular clusters (see Pfahl et al 2009 and references therein) 
the overall rate of such collisions is about 5 — 20 x 10 -10 
per comoving Mpc 3 per year. This corresponds to 10 — 100, 
z < 1 collisions per year. Note that this calculation depends 
in detail on the minimum impact parameter that leads to 
explosive nuclear burning, as studied in detail below. Our 
estimate is significantly lower than a recently submitted pa- 
per (Rosswog et al 2009), primarily due to the assumed core 
radius of globulars clusters. 

Beyond w 50 Mpc, current surveys are not sufficiently 
accurate to distinguish between supernovae in globular clus- 
ters from those that arise from the galaxy field stars in front 
or behind globular clusters (Pfahl et al 2009). Furthermore, 
observations of globular clusters in the nearby SO galaxy 
NGC 7457 have detected what is likely to be a remnant of 
a SNIa (Chomiuk et al 2008). Taken together, these esti- 
mates and observations suggest that white dwarf collisions 
are likely to appear in current and future SNIa samples. As 
double-degenerate supernovae from collisions may not fit the 
standard templates, SN surveys may have to consider pos- 
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Table 1. Simulation impact parameters, 6, d, and total particle 
count, N. 



Scenario 


6[R ] 


d [R w d] 


N [xlOOO] 


1 








800 


2 


0.9 


0.5 


800 


3 


1.7 


0.9 


200 



sible contamination from double- degenerate collisions that 
masquerade as traditional SNIa. 

In this paper we begin to examine collisions between 
two white dwarfs as a pathway for producing SNIa. In ^2] 
we describe our simulations of the impact between two 0.6 
M© white dwarfs, and in ^3] we discuss the resulting hydro- 
dynamics, thermodynamics, and nucleosynthesis. In Q we 
speculate on other mass pairs and impact parameters, and 
discuss the implication of our results. 



2 SNSPH SIMULATIONS 

All of our simulations are conducted with a 3D smooth par- 
ticle hydrodynamics code, SNSPH (Fryer et al 2006). In 
collision scenarios with a non-zero impact parameter, an- 
gular momentum plays a critical role in the final outcome, 
and Lagrangian particle methods generally conserve angular 
momentum better than Eulerian grid methods. We added a 
13 isotope a-chain nuclear reaction network (Timmes 1999; 
Timmes et al 2000; Fryxell et al 2000) and a Helmholtz free 
energy based stellar equation of state (Timmes & Arnett 
1999; Timmes & Swesty 2000) to SNSPH to construct the 
initial white dwarf models and run our collision models. 

Our initial 0.6 Mq white dwarf is composed of equal 
parts 12 C and 16 O and is created using the Weighted Voronoi 
Tessellations method (Diehl & Statler 2006), which arranges 
particles in the configuration corresponding to the lowest en- 
ergy state that is consistent with a given equation of state. 
The resulting white dwarfs have only minor temperature 
variations around lxlO 7 K. Any evolution the white dwarfs 
experience before the collision, then, is due only to their mu- 
tual gravitational interaction. We use a simple 2-body solver 
to calculate the relative velocities of the two white dwarfs 
immediately before a collision for a given impact parameter 
and an initial relative velocity of 10 km s _1 , which corre- 
sponds to the typical virial velocity of a globular cluster. 
The velocities calculated by our 2-body solver then serve as 
the initial conditions for our collision simulations. 

Various angles of incidence for the collision were cho- 
sen. In Table 1, the initial impact parameter, b, is the ver- 
tical separation of the centers of the white dwarfs in terms 
of R at a horizontal separation of oo. Due to gravitational 
focusing, the final impact parameter, d, the vertical separa- 
tion at the moment of collision, is significantly smaller, and 
is given in Table 1 in terms of the radius of a 0.6 Mq white 
dwarf, « 0.01R© . The number of particles N in each model 
is also given in Table 1, i.e., there are N/2 particles per white 
dwarf. 
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Figure 1. Planar slices through the 3D calculation of the evo- 
lution of head-on white dwarf collision scenario 1, given in Table 

1. The features of each panel are described in detail in §3.1. The 
sequence displayed, after the white dwarfs have collided in panel 

2, takes place over a period of ~6 seconds. 

3 RESULTS 
3.1 Scenario 1 

Snapshots of planar slices from the head-on collision case are 
given in Figure [I] This extremely unlikely, but physically in- 
structive scenario was presented in Timmes (2009) and has 
also recently been explored in Rosswog et al (2009). As the 
white dwarfs approach each other, their velocities increase to 
the escape speed, ^4000 km s _1 (first panel). When the two 
white dwarfs first make contact (second panel), shock waves 
develop that attempt to travel outward along the x-axis at 
roughly the sound speed, ~3000 km s _1 (also see Figure [2]). 
The near equality of these two characteristic speeds means 
the shock waves stall as material falls through them (Fig. 
[I] third panel). The region between the two fronts achieves 
a nearly constant temperature (T ~ 10 9 K) and density 
(p rsj 10 6 g cm -3 ) with a non-explosive, mild rate of nu- 
clear burning. This shocked region expands very slowly as 
more material piles into it. In this regime the nuclear energy 
generation rate scales as e w p 2 T 27 . 

When the high density (~ 4x 10 6 g cm -3 ) central core of 
each white dwarf encounters the leading edge of the shock, it 
raises the energy generation rate of enough material above 
the threshold needed to trigger a detonation (see Gamezo 
et al 1999 and references therein). Two curved detonation 
fronts with T ~ 10 10 K form at the interfaces of the shocked 
and unshocked regions and begin to propagate (fourth panel 
in Figure HI which releases enough energy to unbind the 
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Figure 2. The sound speed and infall speed of particles lying on 
the x-axis in scenario 1. This snapshot corresponds to panel 3 in 
Figure^ The material falling into the shocked region is moving 
at or above the sound speed for this region, causing the shock to 
stall. 

merged system (fifth panel). The entire system then under- 
goes homologous expansion (sixth panel), leaving a small 
amount of unburned CO in the central regions surrounded 
by a layer of 56 Ni and other iron- group elements, a layer 
of Si- group elements, and then an outer region of unburned 
CO. The final tally, with 1.2 of CO entering the system, 
is 0.34 M of 56 Ni, 0.33 M of 28 Si, and 0.15 M of unburned 
CO with the remaining mass distributed among other heavy 
elements (also see Figure]!]). Such an event, whose total mass 
is sub-Chandrasekhar, would be a potential candidate for 
either an underluminous SNIa or a new class of transients 
between Novae and SNIa. 



3.2 Scenario 2 

In scenario 2, we use an impact parameter of b = O.9R in 
order to bring the stars into contact with a final impact pa- 
rameter of d = 0.5Rw or d ~ 0.005R© . As in scenario 1, 
the shock front stalls as the inward falling material is trav- 
eling at a comparable speed (first panel of Figure [3]) . Once 
again, when the cores of the white dwarfs enter the shocked 
region, two detonation fronts break out and unbind the sys- 
tem (second panel). However, in this case, the additional 
torque applied to the shocked region as the stars continue 
to move past each other results in off-axis and off-center det- 
onation fronts that occur at earlier times than in scenario 
1 (second and third panels). In this case, 0.30 M of 56 Ni, 
0.37 M of 28 Si are created, leaving 0.25 M of unburned CO 
with the remaining mass distributed among other heavy el- 
ements. 



3.3 Scenario 3 

In scenario 3, we increase the impact parameter to b = 
1.7R , resulting in a pre-collision impact parameter of 
d = 0.9K w d or d ~ O.OO9R . This ensures a purely graz- 
ing incident, which is likely to be the the most frequent 
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Figure 3. Planar slices through the 3D calculation of the evo- 
lution of the d = 0.5K W( i white dwarf collision scenario 2, given 
in Table 1. The features of each panel are described in detail in 
§3.2. 



impact scenario as a wide range of impact parameters will 
result in a similar, grazing collision when tidal effects trans- 
fer angular momentum from passing impactors. Absent the 
violent shock experienced in scenarios 1 and 2, there is negli- 
gible nuclear burning in a purely edge-on impact. Less than 
10 -3 M of 56 Ni are produced by the initial interaction. The 
kinematics of the collision causes both stars to become un- 
bound, forming a rotating disk of white dwarf debris, which 
will eventually collapse into a single compact object as it 
cools. 



3.4 Convergence Studies 

A known weakness of smooth particle hydrodynamic simu- 
lations is resolving shocks with a finite number of particles. 
For this reason, we carried out a convergence study of the 
the d — (scenario 1) and d = 0.5Rw (scenario 2) cases 
in order to place upper limits on the isotope yields. Fig- 
ure^] shows that with increasing particle count, more of the 
CO is converted into 56 Ni and 28 Si, while the abundance of 
32 S remains relatively constant. While we did not achieve 
an absolute convergence in our simulations, we can reason- 
ably extrapolate an upper limit of O.4M of 56 Ni for the 
head-on merger scenario and 0.35 M in the d = 0.5R w d sce- 
nario, where 1.2 M of CO enters the system. In low impact- 
parameter collision scenarios involving larger masses for the 
constituent white dwarfs, it is very likely that this yield 
would increase. 

Rosswog et al. (2009) find 0.32 M of 56 Ni for the head- 
on collision of two 0.6 M white dwarfs with an SPH resolu- 
tion of 2xl0 5 particles. This is consistent with our highest 
resolution SNSPH model with 8xl0 5 particles, but is in- 
consistent with our lower resolution SNSPH models. Ross- 
wog et al. (2009) do not present a convergence study of their 
SNSPH models, and find significantly less 56 Ni is synthesized 
when they use the finite difference code FLASH (Fryxell et 
al. 2009) . We are currently using FLASH to calculate the nu- 
cleosynthesis for the three impact parameter cases studied 
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Figure 4. Isotope yields vs. total particle count for scenarios 
1 and 2. The solid lines are the isotope yields of the head-on 
simulation, while the dashed lines are those of the d = 0.5H w d 
simulation. 



in this paper to investigate this apparent discrepency and to 
validate (or falsify) the 56 Ni yields from our SNSPH models 
(Hawley et al. 2010). 



4 DISCUSSION 

Direct collisions between white dwarfs offer an unex- 
plored mechanism for SNIa production within the double- 
degenerate family of progenitor models. Our simulations 
suggest that at low impact parameters, the collision between 
the most common white dwarf masses (0.6 ) can result 
in explosive nuclear burning, even though the total mass 
involved is below the Chandrasekhar limit, and producing 
enough 56 Ni to result in a dim SNIa or be an example of a 
new class of transients between Novae and SNIa. It is likely 
that collisions between more massive white dwarfs with im- 
pact parameters less than a white dwarf radius will result 
in brighter events. Although our simulations are unable to 
resolve the exact physics that trigger detonations in an open 
environment, these details might not be critical for the over- 
all nucleosynthesis in collision scenarios. 

Rosswog et al. (2009) have reported that a low resolu- 
tion 6 = collision of two 0.9 M white dwarfs can produce 
enough 56 Ni to resemble the nuclear yields of a typical SNIa. 
As we've shown above for two 0.6 Mq white dwarfs, the nu- 
clear yields are sensitive to the impact parameter. Grazing 
incident collisions between more massive white dwarfs are 
likely to lead to similar configurations. The ultimate fate of 
these systems will, like their binary merger cousins, depend 
on the total mass, thermal and neutrono cooling rates, angu- 
lar momentum transport rate, mass accretion rate onto the 
hot central object, and any residual nuclear burning (Yoon 
et al. 2007) . Further studies are needed to say which white 
dwarf collisions lead to double-degenerate la supernovae or 
perhaps a new class of transients between Novae and SNIa. 
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